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Abstract 


The use of superconductors in the field windings of 
large Ssyncnronous machines displaces iron as a major 
element in the production of tne magnetic field and 
promises smaller, less costly machines for a given rating. 
The possibility of replacing the iron in a transformer 
bank with a rotating superconducting field is discussed. 
By placing the primary and secondary three phase windings 
in separate armature windings built in epoxy and fiber- 
glass concentric shells around the superconducting, free- 
Spinning field and allowing the angular alignment of one 
of the armatures to be varied, a unique transformer is 
devised. Unigue features of the synchronous transformer 
are: the ability to separately control real and reactive 
power flow, thus providing power factor correction while 
Simultaneously acting as a transformer, and the ability to 
transform rated real power between two power systems that 
are out of time phase, with zero (or at a specified level 
within a range) reactive power flow. 
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EN ERODUC ELON 


Synchronous machines are commonplace in the modern world 
and the theory of their operation is well developed. The 
application of superconducting technology to the field wind- 
ings of synchronous machines is a new and growing area of 
research [1]. It 1S expected that superconducting machines 
With ratings of 1000 KVA and greater may have as much as an 
order of magnitude reduction in weight and size over a con- 
ventional machine of the same rating with only about, 0.1 
per cent of machine rating required to power the refrigeration 


mt OL the superconductor [2]. 


Recent work at M.I.T. anticipates that it may be pos- 
Sible to manufacture and ship synchronous generators with 
superconducting field windings with ratings as high as 
10,000 MVA and that construction and operation of these 
machines may be more economical than current conventional 


machines [3]. 


These economic advantages are brought about in part by 
the very high flux densities with minimal power dissipation 
achievable with superconductors, which makes it possible to 
eliminate the magnetic circuit except possibly as an external 
shield. Thus an expensive magnetic-iron rotor forging and 
stator teeth may be eliminated, which in turn reduces insula- 
tion requirements (or allows an increase in voltage rating) 


and at the same time increases the amount of space that may 





be filled with armature conductors. Elimination of most of 
the iron from the machine reduces per unit reactances which 


may improve stability. 


References [2] and [4] derive expressions for electrical 
parameters of synchronous machines without magnetic shielding. 
References [5] and [6] extend these equations to machines 


“men an exterior iron shield, 


In a superconducting alternator the superconducting 
field winding displaces iron as a major element in the pro- 
eeteton Of the magnetic field. It might be asked if this 
can be done for other devices, specifically, is it feasible 
to replace the iron ina transformer bank with a rotating 
SM@@ereonducting field? Such a device, hereafter called a 
synchronous transformer, would have two armature windings 
(one for each of two sets of three phase terminals) of the 
form contemplated for a superconducting alternator. It would 
also have a superconducting field winding, which is free 


Spinning as does the field winding of a synchronous condenser. 


Electric energy is transferred from one armature to the 
other by the interaction of the two armature flux waves with 
the field winding flux wave: one interaction being similar 


to that of a motor, the other like a generator. 


As it will be shown, this proposed device has similari- 
ties with conventional transformers, but it has important 


differences too. For example, if the angular alignment of 





one of the two armatures 1S made variable then the device 
has two "handles" on power flow: armature position and main 
field current. This contrasts with the availability of 

Smiy One “handle (tap changing) with conventional trans- 
formers. This thesis will show that these two control 
variables will enable the synchronous transformer to control 
Separately both real and reactive power flow. Thus it can 
act as a power factor correction device while simultaneously 
acting aS a transformer. Additionally, it can transform 
rated real power between two power system buses that are out 
of time phase with no reactive power flow or within a range 
of desired reactive power levels. 

Such a device could be constructed by embedding the 
armature windings in epoxy and fiberglass to a fiberglass 
(or other non-magnetic material) pipe for structural strength 
[7,6]. Two such armature windings would be required. The 
axis of the outer armature would be adjustable so that it 
could be rotated about the longitudinal axis of the rotor 
(and thus also about the longitudinal axis of the inner 
armature) by an external control, but it would not be free 
to spin. The returning flux path and shielding would be 
provided by a concentric iron shell external to the two 
armature shells. Figure 1 illustrates the physical configu- 
ration of the synchronous transformer. Figure 2 illustrates 


eae electrical configuration. 
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Electrical Configuration 


Armature 2 is Shown Rotated by an Angle ¢ from Armature l 
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There is a growing interest in the use of supercon- 
ducting machinery for marine propulsion. In 1966 Sibley 
et al [8] estimated some potential advantages of super- 
conducting electrical machinery over conventional propulsion 
plants. In a later work Laine [9] considered the possibil- 
ity of using a superconducting transmission link in a 
destroyer-type ship. Matjasko [10] investigated a control 
system for a destroyer propelled by superconducting electric 
machinery. Greene [ll] examined transmission systems 
including a cycloconverter scheme between a synchronous 
generator and an 18,000 SHP (13.5 MW) synchronous motor 
iamveng a propeller. 

While in its present configuration it would seem that 
the synchronous transformer would perhaps be more useful to 
the electric power industry than to marine propulsion, it ‘is, 
however, Similar in configuration to a class of dual arma- 
ture Superconducting machines which may have great potential 
as an electric propulsion ship system. This class of AC 
machines may have the control characteristics of a DC machine 
system but with the additional advantage that a solid state 
Switching circuit will replace the bulky cycloconverter. 

This class of dual armature machines will not be discussed 
further here as they are the proprietary information of the 


Smayogenic Engineering Laboratory at M.I.T. 








Note: Unless otherwise noted, 


Table 1 


GLOSSARY OF TERMS 


Ae 


the following subscripts 


will serve as modifiers to the symbols. 


Susscripts 

a phase a, armature 1 

b phase b, armature 1 

c phase c, armature l 

a phase a, armature 2 

B phase 8, armature 2 

Y phase y, armature 2 

d direct axis, armature l 

q quadrature axis, armature 
6 direct axis, armature 2 

2B quadrature axis, armature 
ie field 

O zero sequence, arinar Ure 
02 zero sequence, armature 2 
Py phase variables, armature 
Pp. pnase variables, armature 
Ry rotor variables, armature 
Ro rotor variables, armature 
Ly armature 1 to field 
L, armature 2 to field 
ty equivalent at the terminals, 
e equivalent at the terminals, 


armature l 


armature 2 








Symbols 


Glossary of Terms ae 


(continued) 


GUBLENE 
voltage 

i lbievealigye; te dbi6),< 
inductance 


argulareDOSsTELOn Of pene rOLOr FEOM Ene axis OF 
phase a, armature 1, positive in direction of 
rotation 


angular position of the axis of phase a, 
armature 2, from the axis of phase a, arma- 
Pueewley POSl tt VeuLn the GE Lectl1on-OF rotation 


maximum value of mutual inductance between 
phase a, armature 1, and phase a, armature 2 
ona Loo at ¢ = Q) 


maximum value of mutual inductance between 
phase a, armature 1, and the field 

(dager 5 Lie at 6 = Q) 

maximum value of mutual inductance between 
phase a, armature 2, and the field 

(ce or Lig at 6 = 6) 

Spatial angle from the rotor axis to the arma- 
ture flux wave of armature 1, positive in the 
Q@ipeceion Of retaeloen. 


Spatial angle from the armature flux wave of 
armature 1 to the armature flux wave of 
armature 2, positive in the direction of 
GOudelon 

terminal voltage of armature 1 in per unit 
terminal voltage of armature 2 in per unit 


per unit reactance 


air gap voltage 





uns, 


RESULTS OF A STATIC ANALYSIS 


For any winding of a synchronous machine the voltage at 
the terminals 1s equal to the product of the winding resis- 
tance and current through the winding plus the time rate of 
change of the winding flux. For any winding the flux is the 
Summation of the self inductance times the current through 
the winding and the product of mutual inductances with 
other windings (which may depend on rotor angular position) 
and the currents through those windings. 


thus : 


Vo=rii + (1) 


(2) 


Losses in large synchronous machines are less than one 
per cent of machine rating and the assumption of a lossless 
model does not have much effect on the Sea eae analysis. 
Therefore one can assume that winding resistances are 


neglible yielding: 


dn 


—- —m 
Vee 8a -) 


The machine flux current relationships may be expressed 


by : 
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Ms a Liab Lac HA aie Lay Lif 1 
Ap ety Sige lays, | “Sey be oye 1h, 
Me Lic Linc Le rete Lg Ley Loe Le 
ol See ee Lig Lay Lie iy (4) 
rg Soper ce age =) ey) | Let 18 
A a a a NE PR OR i, 
A bee “oe Yee Use Yee “ye Ye te 


Because of fixed geometry the following relations apply: 


Lew (5) 
L, = lg = L,, (6) 
Lab i Lac 7 Lic ) 
Mog ~ Bey ~ May o 


It is important to note that the geometry between arma- 
tures 1S not fixed since armature 2 may be rotated and since 


phase axes are 120° apart. 


t 50 7 Lg - Bay oF 


It 1s assumed that all rotor to armature mutual 
inductances vary with the cosine of rotor position,with the 
maximum values of Li- and Lie when the rotor axis aligns with 
the respective phase axes of armature 1 and armature 2. This 


is the same as neglecting higher order space harmonics in the 
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Pemeer air gap flux and SESmMAS a non-Salient rotor. 

It is further assumed that armature 2 to armature 1 mutual 
inductances vary as the cosine of the relative angular posi- 
tion of the respective phase axes with the maximum value of 
Low when the armature displacement angle, ¢, equals zero. 
This value Lia (maximum) will be denoted by M. This assump- 
tion is equivalent to neglecting higher order space harmonics 
in the air gap flux between armatures. Higher order space 
harmonics die away more rapidly than the fundamental. Since 
this is an air gap machine, disregarding higher order space 
harmonics in the air gapsS will introduce only a small error 


[6]. Substitution into equation (4) yields equation (10). 


E@Uation (3) gives: 


a a 
V, dp 
ve do 

2 d 

Yo | = at | *o aa 
y r 

g g 
V r 

uf Y 
Ve he 


At this point it 1S convenient to transform equation 
tei} into the rot2zting coordinate system of the rotor. The 


Subscripts d and q correspond to quantities of stator l, 
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petating With the direct and quadrature axis of the field. 
The subscript, o, corresponds to zero sequence quantities of 
stator 1, that is in phase in all three windings of stator 1. 
The corresponding subscripts for stator 2 are 6, &, and °. 
The details of this transformation are given in Appendix A. 

The result of transforming to rotor coordinates yields 
eamectons (12) and (13). 

By comparison of equations (12) and (13) the following 


relationships can be recognized: 


Ci Wg (14) 
Ve ~ = (A) + WAS Ge) 
Ve = & (Ag) - wre (16) 
V, = a (Ap) + wr (17) 
Ve = a¢ Ag (18) 


It can be seen that zero sequence variables can be 
elemimated from equations (12) and (13) since they are com- 


pletely uncoupled from other variables. 


For ease in notation the following definitions are 


introduced: 


(ee), 


(20) 


Thus equation (13) becomes: 











d E 

rg 0 

ds = = M cosd 
he 5 M sind 
re - af 


8 
0 aa cos¢ 
L 2 M sind 
E Ph 
3 
5 M sind oe 
3 
5 M cos¢d 0 
/3 
ols. 


Ko] Go 





= M sind 7s Ly 1, 
M cos¢ 0 - Ph 
/3 Uo) lig 

0 i; 

Be a, 

ZL} 


Now that the basic relations governing the synchronous 


transformer have been developed, the operation of the trans- 


former between two buses will be looked at. 
vention will be observed throughout, 


input to the machine from either system (bus) 


that is, 


Motor sign con- 
real power 


is defined as 


positive and lagging current determines the positive direc- 


tion for reactive power. 


Since the flux-current relationships have been trans- 


formed into rotor d-q coordinates, 


the following definitions 


of armature equivalent flux vectors are introduced: 


eimilanrly, 


—_ 
hyo VX 5 + he 


(22) 


(23) 


the equivalent armature voltage vectors are 


likewise defined: 


oe el 
0 ee 
a5 ay 


(24) 


(25) 








23 
The angle o is defined as the Spatial angle from the 
rotor axis to the armature flux of armature 1, positive in 
the direction of rotation. The angle 6. is defined as the 
spatial angle from the armature flux wave of armature 1 to 
the armature flux wave of armature 2, positive in the direc- 
tion of rotation. Figure 3 illustrates these sign conven- 


EOS » 


Pakcetndeively, since the voltage vectors are a ninety 
degree rotation of the armature flux wave vectors, these same 
definitions can be viewed as applying egually between the air 
gap voltage, terminal voltage of armature 1 and terminal 
voltage of armature 2. Figure 4 illustrates the sign conven- 


tions of angles 6 and oe in terms of the voltage vectors. 


mae tollLowing definitions follow directly: 


Veuve noe (26) 
Ne = Vi, cosé. (27) 
Vs aah rie sin(o, + 6.) f2c)) 
Ve eS cos(6. + S50 “G2o)) 


If the restriction to strict steady state operation is 
assumed, then time derivative terms become zero. Equations 


(14) through (17) become: 
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Figure 3 


Illustration of Sign Conventions 
for Spatial Angles 6 and oe 
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Figure 4 


Hmblustration of Sign Conventions for Spatial Angles 
s and 0, in Terms of Voltage Vectors 





tO 
Or 


Vara Wr (30) 
Ma = wha (31) 
Ve WH (2:2) 
Vy = Wr. (338) 


Substituting of equations (21) through (29) into 


equations (30) through (33) yields: 


oa ; S ee 5 , 
Vid cosé. = ee + 5 M sinoi, + 5 M cos$i;) (34) 
V Simo. = w{L 1, + 2 M cos¢l1, - 2 M sindi, + /3 ie eas oe 
sgl a Ed 2 6 2 2 V2 ae 
V.5 cos(6_ + 6.) = -w(- 2M sindi, + 2M cossi_ +L i-) (36) 
12 a O Zz d 2 q oaG 
V Sato +6) = ie M coso1, + 3 M sindi. +Li, + /3 de a 
tea ayy 0 2 d 2 q [ag De 
(37) 
Since the rotor is free spinning there will be no 
electrical-mechanical power transfer from the rotor. If we 


ignore all losses (friction, windage, copper losses, etc.) 
then the net power of the synchronous transformer will be 


zero. 


0 = V 3 Vie Ved cee Vea (33) 
qq 


ata E 


This is not to say that the device cannot absorb power 


from one bus under the lossless assumption, but rather 








equation (38) tells us that real power into either 27 
armature must go out the other armature. 

Substatution Of equations ((14) through (17) with the 
Steady state asumption d/dt = 0 and equation (21) into 
equation (38) yields: 


0 = AEs 2s - Whois: (39) 


Equations (34) through (37) and equation (39) are the 
governing equations for the synchronous transformer operating 


between two buses. Independent variables are i >, V 


ff? ul 
and V5: Dependent variables are Las loss dat Ley our and So: 


Thus we have five equations and six unknowns. The method of 
solution used is to assign values of Si for a eee 
Wide range and at each value of - solve the resulting equa- 
tions whicn then have only five unknowns. However, before 
proceeding to the solutions of these equations it is con- 
venient to put them into the per unit system. Details of 


this are in Appendix B. 


The governing equations for the synchronous transformer 


Operating between two buses in per unit are: 


Vy sino. = Behe + x1 6 Sing + ete cos$¢ (40) 

Vy cose. = Xolg + X16 cos®o - Xie Sing + X ig (41) 

es sin(6 +0.) = - xX ig Sind + ae cos¢ + Xgi- (42) 
V5 cos (6_+5) = x tg cos¢ + Sg Sing + X ois + Xp ols (43) 
(Ol a. come: epic, ae cinerea (44) 








In order to use Gauss-Seidel computer aided iterative 


solution these equations were rewritten as follows: 





ae years a Xt s Sai a Xmte cos¢ | ae 
- = Sin Vy 
. Vy cosé. - X)15 cos¢ + Xnte Sind - Xr te 
Lg Sg eas (46) 
se 
d 
fe. V5 sin(6 +0.) + x tg Sind - ane! cos¢ 
L_ = —--oO oo (47) 
e me 
6 
em . V5 cos (+6) - x ig cosd - Xnig Sino - Xp ole (48) 
5 x 
6 
aaa a 
i,=- at (49) 
a sap 


Appendix (C) discussed the Gauss-Seidel iterative solu- 
tion. Reference 12 discusses this and other applicable 


methods in great detail. 


Reasonable values of per-unit inductances for a syn- 
chronous machine with a superconducting field winding were 
chosen [6, 13]. The assigned inductances are listed in 


Table 2. 








Table 2 


Bom UNt TT REACTANCES USED IN STATIC ANALYSIS 


Armature 1 Synchronous Reactance Xa = Op au. 
Bimatuke s2 Synchronous Reactance Xo = Oe > pre Ue 
Mutual Reactance Between Armatures x, = 0.4 p.u. 
Armature 1 to Field Resistance X11 = POS rule 
Armature 2 to Field Reactance X1> = LO peu. 


A computer program was prepared to solve simultaneous 
equations (45) through (49) by the Gauss-Seidel eee 
method for a range of armature displacement angles, $, from 
-90° to +90° and a range of spatial angles between armature 
voltage vectors, Bor from 0° to 90° and for different levels 
of field current. In all cases reactances from Table 2 and 
one per unit terminal voltages were used. For each level of 
field current and each choice of angles $6 and - the computer 
calculated and printed values of direct and quadratic axis 
currents and voltages for both armatures, real and reactive 
power for both armatures, and the spatial angle between the 
air gap voltage vector and the voltage vector of armature l. 
Tnis produced a wealth of data from which Figures 5 through 


10 were extracted. 


It is important to note here that while the angles 5, 
and 6, are spatial angles relating the vslitage vectors to the 
rotor, the angle ¢ 1S a physical angle by which the armatures 


are displaced.. Thus the time phase angle between the power 





50 
systems attached to the armatures iS a combination of their 
physical and electrical angular difference, - Oo Nee 
we are not interested in internal angles to characterize the 


device this power system time phase angle will be the sig- 


nificant angle used to characterize the device. 


Since we are primarily interested in the synchronous 
transformer to transfer power, the solutions yielding real 
power equal to one per unit were collated and graphed in 
Figures 5 through 7. Figures 8 through 10 were constructed 
from Figures 6 and 7. Figure 8 was obtained by cross plot- 
ting the armature displacement angle, $, from Figure 6 
corresponding to the zero crosSings of reactive power flow 
Peomerrgqure 7. Thus Figure 8 gives the setting points for 
~eiweiee displacement angle, $¢, and field current, I-, for a 
range of power system time phase differences in which one per 
unit real power can be transferred with no reactive power. 
Figures 9 and 10 yield similar results for reactive power 


Pevmequal to plus and minus 0.1 per unit. 
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Figure 5 


Curves of armature cisplacement angle, ¢, versus 
Spatial angle between armature voltace vectors for 
one per unit real power transfer and various levels 
Oberle lorcurrent, I. 
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Figure 6 


Curves of armature displacement angle, ¢, versus 
power system time phase angle difference, 6, + $, 
for one per unit real power transfer and various 
levels os field current, I, 
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Curves of reactive power, Q, versus power system 
time phase angle difference, 6. + ¢, for one per 
unit real power transfer and various levels of 
ficla Current, I, 
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Figure 8 


Armature displacement angle, $6, and field current, 
iene veseus pewer ssyotem tine piace vangle diiiers 
ence, 6, + ¢, for one per unit real power transfer 
and reactive power egual to zero. 
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Figure 9 


Armature displacement angle, 6, and field current, 
We eversus POWer System, cime Phase angie Critter — 
eens os. | 7, LOG One per Unie meal power _ranster 
and reactive power equal to + 0.1 per unit. 
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Armature displacement angle, $¢, and field current, 
I-, versus power system time phase angle differ- 
ence, 55 + ¢, for one per unit real power transfer 
and reactive power equal to -0.1 per unit. 
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DiSCUooLON Or eeoouUiseS 


The important result of this thesis is that the 
synchronous transformer can be characterized by a set of 
curves in which the armature displacement angle and field 
current can be determined for a given power system time 
phase angle difference for operation at real power equal to 
one per unit and at various levels of reactive power. In an 
ordinary transformer there is always reactive power flow due 
to leakage reactance, but in the synchronous transformer 
Settings may readily be chosen for unity power factor power 
transfer. Additionally, the synchronous transformer can 
Operate at unity power factor wien transferring power between 
two systems that are out of time phase or it can simul- 
taneously act as a transformer and a synchronous condenser. 
between two buses that may be in or out of time phase. Thus 
it is seen that a synchronous transformer in addition to 
being able to act as a transformer and as a synchronous 


condenser, can also act as a phase shifter. 





CONCLUSIONS AND RECOMMENDATIONS 


@enclusions 


The results of this static analysis have been encourag- 
ing. It has been shown that the analytic model of equations 
(40) through (44) yield a set of curves in which the armature 
displacement angle and field current can be determined for a 
given power system time phase angle difference for operation 
at real power equal to one per unit and at various levels 
of reactive power. The synchronous transformer Sein eraiag eau 


real power from one bus to another whether its time phase 


fmmehind, equal to, or ahead of that of the other bus. 


Recommendations 


It is recommended that further investigations of the 
Syiienronous transformer be pursued. Further efforts should 
include a dynamic analysis with faults and transients and 
an economic analysis to determine tne best level of its 
usage within the power systems and its general economic 


feasibility. 
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APPRENDI te 


TRANSFORMATION TO ROTOR COORDINATES 





TRANSFORMATION TO ROTOR COORDINATES eos 


The machine voltage current relationships (equations 


(10) and (11) may be expressed by equation (A-1) (page 43). 


Define a pair of reciprocal transforms relating phase 


and rotor corrdinates, where u represents a current or a 


voltage. These transforms are similar to those of Lewis [14]. 
Us u. 
te = 4 UL, (A-2) 
Yow LY 
u. Us 
2 eel _ 
Uy, = T U4 (A-3) 
bare lige 
~ cos@ .cos(@ - 27/3) cos(6@ + 21/3) 
T = /% Seine —sin(e no 1 7s) on cme 72) (A-4) 
1/2 az V2 
cos8@ -siné 1 
T = /% Gos (8 =| 20/3). estan 2 el (A-5) 
e6s(6 -4+°21/3))) =san (cee 2, 3) ee 


By partitioning equation (A-1) as shown we can reduce 


the notation as shown below and treat the partitions 


separately. 
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Note that these submatrices are not symmetric, i.e., 


A 7 = iS - The off diagonal elements are transposes of their 


Counterparts. 


The transformed version of equation (A-6) is: 





d -1 d -1 d a 
Vel Sa eee ae Bae oe | | ea 
: a = d =a 
Veo i: eae ate oa" Bae A: tr? Jaa) 
d -l1 a -l a ays 
Ve ES ‘ea at fot =a | (tara 


Carrying out the above matrix operations yields equa- 


meen (12) of the text. 


If the flux-current relationships (equation (10)) are 
partitioned in the same manner as equation (A-1), they 


become: 


oa a A ao Bat igatk 
ND Ure lacie eo: Aye| |*p2 ESS 
Mech ita “ta Aeaps fre 


The transformed version of equation (A-8) is: 
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Carrying out this matrix operation 


(13) of thetext. 


(A-9) 


yields equation 





APPENDIX B 


THE PER-UNIT EQUATIONS 





THE? PER-UNITT VEOUATIONS 


There are various techniques for changing equations into 
per-unit equations, all of them related. The method that I 
Shall employ parallels that of the notes for Power Systems 
Peemmeering Il, course 6.552 at M.I.T., by Prof. G. L. 
Wilson. 


The following base quantities are specified: 


Ax1s Base Voltage Base Currene 
d V if 
O O 
q Vo ve 
Noy wa 
° Yo N. to 
a Ql 
Ny Na 
é Vo N_ sO Ne 
a QQ. 
: cae nea 
ON Zee 
a fe 


The subscripted N's in the above base quantities 


represent turns ratios. 


tie equations tO be put inte per un? cease the sgoveruine 
equations of the synchronous transformer operating between 
two buses, equations (34) through (37) and (39). These are 


listed below: 





3 3 ; 
V, cos on Neth + 5M sing 1, + 5M cos@¢ 1 


¢) 


o | 3 oe Ale [a 
Vi sin = = w (Lots + 5M cos$é i 5M Sind if + WS Lji-g) 
(Ba2) 
ee ees ; 3 7 
V5 cos (0 .+5)) = —w ( aM sind 1a + 5M cos¢é te s bere) (B-3) 
V. sin(é6_+6_) = w (2M COS ae 34 Sind i. +Li, + /3 pak ed 
2 VS. 2 d 2 gq ec 6 y awe ae 
(B-4) 
0 = cole a oe (B-5) 


However, we recall that equations (B-1) through (B-4) 
came from equations (14) through (17), which in steady state 


per-unit become: 


Wal puus) WGN cUpaut ae 
ye) eee o OF aes ea) 
Oia itas Bete = re We) os (B-8) 
eee a eee (B-9) 


Thus we see that the quantities which we must find in 
per-unit are the fluxes. The flux equation is equation (21), 


which is broken down into individual flux equations below. 
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TS 3 N° 
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Now identify the following quantities: 
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Now substituting equation (B-22) through (B-26) into the 
foregoing flux equations and dropping the parentheses and the 
Subscript "p.u." since all quantities are now in per-unit 


yields: 





ay IL 


h 3 = Xalg + x1» cos¢ - Xie Sino + Xr le (e270 
a = as + x 1s Sind + Xnte cos¢ (BaZa 
ds = xt cos¢ + eae Sind + Xai. + Xp ots (B-29) 
ry =-X 1, Sind + ee cosd + Xie (Ba30) 


Substituting these results into equations (B-6) through 
(B-9) and noting that for steady state synchronous operation 


W/W mee D.sUs yrelds: 


Ve = Soe + x1 3 Sind + te cos $6) (B-31) 
Ve = Xals + xX ts cos¢ - Xnte Sind + Xr te (B-32) 
V5 = Xt Sind - ae cos¢ - Xgl; (B-33) 
ie = x14 cos¢é + Both Sind + Xl. r Xrol¢ (B-34) 


Similarly, equation (8-5) can be put into per-unit as 


follows: 


U) ; , 
0 = — L, (i_) I + — Lo(1 
Wo I SG sa. 4 6 Wo 2 


Ey pu. tea ole (B-35) 


Multiplication of both terms by the proper constant and 


assuming steady state operation yields: 


Q = cad + Xp ogi; (B-36) 





BauUeat1ons {(B-31) through (B-34) and (B-36) are the 
governing equations for the synchronous transformer in per- 
Umat . 


The following definitions are introduced: 


= Fa Dis ts (B-37) 


we = eo oe (B-38) 


Substitution of the definitions of equations (26) through 
(29), (B-37), and (B-38), directly yields the per-unit equa- 
tions for the synchronous transformer operating between two 


buses. 


V,sind_ - Sal + x1» Sing + te cos¢ pB= 39)) 

V, cose. = X51 =e xX ts cos¢é - ie Sing + Xp ile (B-40) 
Vosin(s_+6) = -x 1, Sind + Sets cosé + Xgl; (B-41) 
V5cos (5 +6.) = X15 cos@é + ae Sing + Xo. S Xp ols (B-42) 
Oy St. See ee (B-43) 





APPENDIX C 


GAUSS-SETDELsMELie» 





GAUSS-SEIDEL METHOD 


The Gauss-Seidel method is a useful technigue for 


solving the following system of simultaneous linear equations: 


Sibi | eq es Se eee, 
eee 0 op) ee re Ome 
ordi = ORES, Cee oe Ga Gnn*n = on ntl (C-1 


iordaer GO Start the algorithm an initial solution 


vector must be assumed and the interative algorithm is: 


*i ntl r ag 
ae (C-2) 
1 aad. ee wee ac J 
pEL = epee 
Eipe 2 


In the iterative algorithm, (C-2), the most recently 
available values of ee are always substituted in on the 
right-hand side. 


The convergence criterion is: 


th 


|x, (Caine iteration) - x, (k iteration)| < « (C=3) 


B= Ves 2 eel 
For computer solution an upper limit on the number of 


iterations should be specified since convergence may not 


occur. 
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